Configuration interaction calculation for the isotope shift in Mg I. 



J. C. Berengu£| and V. A. Dzuba 
School of Physics, University of New South Wales, Sydney 2052, Australia 

V. V. Flambaum 

School of Physics, University of New South Wales, Sydney 2052, Australia and 
Institute for Advanced Study, Einstein drive, Princeton, NJ 08540, USA 

M. G. Kozlov 

Petersburg Nuclear Physics Institute, Gatchina, 188300, Russia and 
Queen's University of Belfast, Belfast, BT71NN, UK 
(Dated: 23 December 2003) 

We present an ab initio method of calculation of isotope shift in atoms with a few valence electrons, 
based on the configuration-interaction calculation of energy. The main motivation for developing the 
method comes from the need to analyze whether differences in isotope abundance in early universe 
can contribute to the observed anomalies in quasar absorption spectra. The current best explanation 
for these anomalies is the assumption that the fine structure constant a was smaller at early epoch. 
We show that we can calculate the isotope shift in magnesium with good accuracy. 

PACS numbers: 31.30.Gs, 31.25.Jf 
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The motivation for this work comes from recent stud- 
ies of quasar absorption spectra that reveal a possible 
change in a since the early universe pj . One of the pos- 
sible major sources of systematic effects in these studies 
is that the isotopic abundance ratios in gas clouds in the 
early universe could be very different to those on Earth. 
A "conspiracy" of several isotopic abundances may pro- 
vide an alternative explanation for the observed variation 
in spectra 0. In order to test this possibility it is nec- 
essary to have accurate values for the isotope shift (IS) 
for the relevant atomic transitions. Experimental data is 
available for only a very few of them; therefore, accurate 
calculations are needed to make the most comprehensive 
analysis possible. 

Previously we have calculated isotope shift in atoms 
with one valence electron 3] . This work represents a first 
step in developing a method for the calculation of isotope 
shift in atoms with more than one valence electron. The 
method used here is similar to our previous calculations 
of the a-dependence of transition frequences for ions with 
several valence electrons Q. It includes Dirac-Fock cal- 
culation of the core and configuration interaction (CI) for 
the valence electrons in combination with the finite-field 
method for the perturbation. Magnesium is one of the 
simplest and well studied two-electron atoms. Because of 
that it is often used as a test ground for different methods 
of atomic calculations. In this paper we show that we can 
calculate the isotope shift of some magnesium transitions 
for which experimental values are available. 



The isotope shifts of atomic transition frequencies 
come from two sources: the finite size of the nu- 
clear charge distribution (the "volume" or "field" shift), 
and the finite mass of the nucleus (see, e.g. |{|). 
The energy shift due to recoil of the nucleus is 
(l/2M)p% = (l/2M)(Sp l ) 2 . Furthermore this "mass 
shift" is traditionally divided into the normal mass shift 
(NMS) and the specific mass shift (SMS). The normal 
mass shift is given by the operator (l/2M)£p 2 , which 
is easily calculated from the transition frequency. The 
SMS operator is [l/M)T,i < j(pi ■ pj) which is difficult to 
evaluate accurately. 

The shift in energy of any transition in an isotope with 
mass number A' with respect to an isotope with mass 
number A can be expressed as 



A' A 



NMS 



l^[±j-j)+F6(rY' A , (1) 



where the normal mass shift constant is 



&NMS 



1823 



(2) 
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and (r 2 ) is the mean square nuclear radius. The value 
1823 refers to the ratio of the atomic mass unit to the 
electron mass. 

In this paper we develop a method for calculating the 
specific mass shift fcsMS for atoms with several valence 
electrons. It is worth noting that in this paper we use 
the convention 8v A ' A = v A — v A . 

Following our previous work on single valence electron 
atoms (Ref. 0) we are looking for an "all order" method 
of calculation. Again we have found that the finite-field 
scaling method is very useful in this respect. The rescaled 
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SMS operator is added to the many-particle Hamiltonian 

Hx = H a + XHsus = H + X ^2 Pi ' Pi • (3) 

i<j 

The eigenvalue problem for the new Hamiltonian is solved 
for various A, and then we recover the specific mass shift 
constant as 

ksms = hm — . (4) 

The operator @ has the same symmetry and structure as 
the initial Hamiltonian Hq (see the Appendix, Ref. 0). 

In this work we restrict ourselves to the frozen-core 
approximation. We first solve the Dirac-Fock equations 
for the core and valence electrons. Then we generate 
a basis set that includes the core and valence orbitals 
and a number of virtual orbitals. Finally we do the full 
configuration interaction calculation. 

The SMS operator for the valence electrons in the 
frozen core approximation can be divided into the core, 
one-electron and two-electron parts: 

H S MS = #SMS + ^SMS + H SMS ( 5 ) 

The first term in Eq. (J5J corresponds to the change of 
the core potential. It accounts for the change of the core 
orbitals when the Dirac-Fock equations are solved for the 
operator H\. The term -ffsMS accounts for the exchange 
interaction of the valence electrons with the core: 

JV coro 

(*l#SMsl fc > = J2(hj\Pi-P2\j,k). (6) 

3=1 

The last term corresponds to the specific mass shift be- 
tween the two valence electrons, p± -p2- 

III. CALCULATION AND RESULTS 

We are using the Dirac-Fock code , which was mod- 
ified for the Hamiltonian ©. The CI calculations are 
made with the help of the modification of the code 
Q . In order to study the role of the valence correlations 
we made three different calculations: 

1. The basic one-configurational calculation for the 
ground state 1 S'o[3s 2 ] and for the 1 ' 3 Pj[3s3p] states. 
All core orbitals and orbitals 3s and 3p are formed 
in the V N ~ 2 approximation (i.e. by solving Dirac- 
Fock equations for the core). 

2. Full two-electron CI for the medium size basis set 
[8sp5d], which includes the orbitals 1— 8si/ 2 , 2—Spj 
and 3 — 5<ij. The 3s, 3p and 3d orbitals are solutions 
of the V N ~ 2 Dirac-Fock potential. The remaining 
virtual orbitals are constructed by multiplying the 
previous orbital of the same partial wave by the 
simple radial function and orthogonalizing with the 
other orbitals 



TABLE I: Calculations of the specific mass shift constants 
fcsMS for Mg I levels (in GHz ■ amu) . Individual contributions 
from Eq. (JKJ are presented, as well as the total. For some 
levels we give medium (M) CI and one-configurational results 
(S) in addition to the large (L) CI ones. 







fcsMS 






Level 


(0) 


(1) 


(2) 


E 


CI 


1 S (3s") 


559 


-883 


131 


-193 


L 


11 


561 


-881 


135 


-186 


M 




857 


-1125 





-268 


S 


3 Si(3s4s) 


422 


-615 


44 


-149 


L 




431 


-624 


52 


-142 


M 


1 S , o(3s4s) 


415 


-615 


21 


-179 


L 


11 


424 


-630 


30 


-177 


M 


1 D 2 (3s3d) 


343 


-616 


-267 


-541 


L 


3 Di(3s3d) 


375 


-561 


41 


-144 


L 


11 


381 


-571 


-10 


-200 


M 


3 D 2 (3s3d) 


375 


-561 


41 


-144 


L 


3 D 3 (3s3d) 


375 


-561 


41 


-144 


L 


3 P °(3s3p) 


428 


-853 


-144 


—570 


L 


3 i?(3a3p) 


428 


-852 


-145 


-569 


L 


3 P 2 °(3s3p) 


428 


-850 


-145 


-567 


L 




431 


-850 


-142 


-561 


M 




759 


-1161 


-266 


-668 


S 




408 


-698 


329 


38 


L 




411 


-700 


341 


52 


M 


11 


946 


-1163 


265 


49 


S 


3 P °(3 S 4p) 


402 


-630 


13 


-215 


L 


3 Pi°(3s4p) 


402 


-629 


13 


-215 


L 


3 P 2 °(3s4p) 


402 


-629 


13 


-214 


L 



3. Full two-electron CI for the basis set [12spc?9/]. 
This basis set is formed by diagonalizing the Dirac- 
Fock operator on the basis set of B-splines and ex- 
cluding orbitals with high energy (for a description 
of this method as applied in atomic physics, see e.g. 
Ref. P3). 

Below we refer to these calculations as small, medium, 
and large. The large calculation is already very close to 
the saturation of the valence CI. Here the difference be- 
tween the theoretical spectrum and experiment is mostly 
caused by the neglect of the core-valence correlations. 
The latter were studied, for example, in Ref. |nj. For 
Mg I the typical corrections to the valence energies and 
transition frequencies were found to be of the order of a 
few percent. 

Table[I]presents the resulting SMS level shift constants, 
^SMS of Eq. in different approximations. The contri- 
butions of individual terms in Eq. (|SJ| are given, as well 
as their sum. It is interesting to note that all of the terms 
are large in comparison to the total SMS. There is a large 
cancellation between contributions within levels, and also 
between different levels. This shows that very high accu- 
racy is required in order to preserve the remaining SMS 
in transitions. 

Comparison of the different approximations shows a 
strong dependence on the size of the basis sets. We see 
that it is very important to saturate the basis as com- 
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pletely as possible. In some cases the SMS changes dras- 
tically even between the medium and the large basis sets. 
In particular, the difference between large and medium 
SMS calculation for the level 3 D 1 (3s3d) is 39%. That is 
mostly due to the /-wave contribution, which is absent 
in the medium basis set. Note that the SMS operator 
can only mix orbitals with Al = 1. That is why the /- 
wave contribution is more important for the levels of the 
configuration 3s3d On the other hand, for the same rea- 
son, the contribution of the higher partial waves to the 
considered levels is suppressed. 

Analysis of Tablc^shows that valence correlations tend 
to decrease the contributions of the first two terms of 
the SMS operator. The third (two-particle) term of the 
SMS operator is generally not screened. On the contrary, 
for some levels the two particle contribution grows with 
the size of the basis set. Note that the final value of 
the two-particle contribution to the ground state SMS 
is of the same order as the other contributions, as it is 
for most other states, while in the one-configurational 
approximation it is zero. 

In Table [H] we compare the results of our calculation 
with experiment for SMS in transitions between 26 Mg 
and 24 Mg. Also presented for comparison are the results 
of Veseth (Ref. |lj). That paper used non-relativistic 
many-body perturbation theory within the algebraic ap- 
proximation to calculate the isotope shift to third order 
for some transitions. 

We have also applied the finite-field scaling method to 
calculate the field shift in Mg. By following the defini- 
tions for field shift given in Ref. Q, and using the same 
approximations for the CI calculation, we found that it 
was less than 2% of the normal mass shift for all relevant 
transitions. The field shift is smaller than omissions in 
the mass shift calculation, notably the core-valence cor- 
relations; thus for simplicity we have neglected the field 
shift from our analysis in this paper. 

Core-valence correlations have been studied for the 
one-electron atoms in Ref. 0, ll3T | and shown to be quite 
noticeable. They can explain the difference between our 
calculations and the experiment in TablelTTI Core- valence 
correlations are usually more important for the ground 
state than for excited states. That may be the reason 
why the largest discrepancy with the experiment is for 



the transitions from the ground state. 



IV. CONCLUSION 

We have presented a method for the calculation of 
the isotope-shift in many-electron atoms using the CI 
for the valence electrons in combination with the finite- 
field method, and tested the method in magnesium. The 
agreement was found to be quite good for all transitions. 
Even for the transitions from the ground state 1 Sq to the 
J = 1 levels of the configuration 3s3p, where the error is 
largest, it constitutes about 20% of the total IS. In par- 
ticular, for the purposes of resolvingsystematic errors in 
the search for a- variation (Ref. [l|,|2|), such accuracy is 
high enough. 

Further work on magnesium could include core- 
correlations, using the extensions to CI outlined in 
Ref. 0- We have decided not to do this, however, be- 
cause we wanted a general method for calculating IS in 
many electron atoms. The method of including core- 
valence correlations in the valence CI with the help of the 
effective Hamiltonian has proven to be very effective for 
atoms with two or three valence electrons, but becomes 
less reliable for atoms with more than three valence elec- 
trons. Unfortunately, most of the ions of astrophysical 
interest have many electrons in the open shells. For such 
ions valence correlations are the most important ones and 
we plan to use this technique to calculate isotope shift for 
the transitions that were used to detect variation of a. 
That will provide stringent limits on the size of the sys- 
tematic error due to variation in isotope abundance. 
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TABLE II: Comparison with experiment of the SMS for several transitions (in MHz) between 26 Mg and 24 Mg. Also presented 
are the results of Ref. j "1 21 ] for a theoretical comparison. We have assumed that the field shift is negligible. 









A 


IS (expt.) 


NMS 




SMS 


Transition 


(A) 






Expt. 


Present Ref. [12] 


i S'o(3^) 
1 S (3s 2 ) 


— > 


*Pi°(3s3p) 


4572 


2683(0) a 
1412(21)' 


1153 


1530 


1205 1378 


-» 


'pr^p) 


2853 


1848 
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-740 






11 




1390(31) c 




-458 




3 PS(3s3p) 




3 Si(3sAs) 


5169 


-396(6) d 


1020 


-1416 


-1349 


3 Pi°(3s3p) 
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5174 
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-390(7) d 


1019 
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3 P 2 °(3s3p) 
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3 Po(3p 2 ) 
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1810(80) E 
60(15)' 
61(3) ' 
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-1316 
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3 P?(3s3p) 
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3 P 2 °(3s3p) 




3 D lj2 , 3 (3s3d) 3839 
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-1315 
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3 P?(3s3p) 
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3094 


420(20) e 
2107(15) d 


1704 


-1284 


-1375 


'pr^p) 




1 D 2 {3s4d) 


5530 


953 


1154 


1224 
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